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Q \ Abstract 

, We make a systematic investigation of the generic properties of static, spherically symmetric, asymptoti- 

cally flat solutions to the field equations describing gravity minimally coupled to a nonlinear self-gravitating 
O , real scalar field. Seven corollaries and a theorem on selection criteria for two- and one-parametric solu- 

tions are proven and conditions for obtaining particle-like solution, black holes or naked singularities are 
5-H , derived. A series of exact solutions in closed forms describing different black holes, naked singularities and 

particle-like solutions are provided. 
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Due to the accelerated expansion of the Universe, scalar field dark matter models emerged to account for the 
astronomic observations (see [1] for a review). The models describe dark matter as a scalar. These models 
^ ■ are special issues of more general theories known as scalar-tensor theories. In contrast with pure gravity, four- 
dimensional scalar-tensor gravity admits a variety of static, spherically symmetric, neutral or charged configu- 
rations including Schwarzschild black holes (BH), phantom BHs, naked singularities (NS), particle-like (P-L) 
^ ■ solutions and wormholes [2-8]. 

We consider the four-dimensional action of general relativity coupled to a self-gravitating real scalar field 
with positive kinetic energy (Einstein-scalar gravity) 



S = Jd\^[R + (j),^^'^-2Vm, 



(1) 



where R is the scalar curvature and V ((^) is a potential for (^. The field equations derived from ([TJ are [5] 



fV^V^0+y' = O, R%-^5^yR + T%=0, (2) 



(where ' = d/dp), with T^^y being the conserved energy-momentum tensor of (p 



For static, spherically symmetric configurations, the metric can be put in the form 

dp 2 

m 



ds^ =A{p)At'^ -^^-?{p)AQ} . (4) 
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Figure 1: Typical behavior of r(p). c is a scaling constant and location of the center of configuration. 

We fix the behavior of r(p) at spatial infinity such that r ~ p as p ^ oo (this is no loss of generality). With the 
ansatz (01), the first equation in ^ reduces to Eq (O and the second group of equations in ^ lead to Eqs (O 
to ^ [5] 

rV = f(Ar20')', (5) 
V = -{A'r^y/{2r^), (6) 
0'2 = _2r7r, (7) 
A{r^)"-r'-A" = 2, (8) 

Only three equations of this set are independent; when the last three equations, (O to ([8]l, are satisfied the first 
one is automatically satisfied. 

Theorem 1 of Ref. [5] and its generalization to {d > 4)-dimension [6] rules out wormhole, horn and flux 
tube type solutions to Eqs ^ to ([Hi. A consequence of that, the range of r is from r = at the center to oo; 
correspondingly, the range of p is from p = c at the center of configuration (r(c) = 0) to oo. Since <p'^ > (real 
scalar field), ^ leads to r" < 0. Hence, the function r(p) is concave down on its range. Fig.[T]shows a generic 
typical behavior of r(p ) where c is the location of the center. In this paper c is treated as a(n arbitrary) scaling 
constant (which can be set to 1) to which we compare the masses of the solutions and the whole theory is scale 
invariant, c is also a measure of the strength of the field. We also make it such that the value of p at r = is 
precisely c; this is no loss of generality as justified in Appendix A. 

For the case of four-dimensional field equations ^ to ([D, Bronnikov [5] established some restrictive results, 
which have been latter generalized to J-dimensional field equations [6]. We are concerned with the following 
conclusions [5,6]. 

• All horizons are simple. 

'Setting c = 1 is the same as the rescahng p/c ^ p, m/c ^ m, A A, (j) ^ (j), (?V V. 
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• All possible global structures for (j) ^ constant and for (j) = constant are the same; 

• Possible non-singular solutions have either Minkowski/ Ads or de Sitter structures. P-L asymptotically 
flat solutions with mass m > and V having both signs are not excluded; however, if V > no such P-L 
solutions (m > 0) exist. 

• For all functions V{p), Eqs (l5]l to ([8]l do not admit static, spherically symmetric regular BH and wormhole 
solutions with any asymptotic. 

• For V >0, Eqs Q to ([H) do not admit any hairy BH solutions (no scalar hair theorems [9]). 

An asymptotically flat BH solution to Eqs ([5]l to ([8]) is characterized by the presence of event horizons. Since 
all horizons are simple, there is a change in the sign of A(p) at each event horizon. Thus for a BH, A(p) has 
both signs in the domain c < p < oo. For an asymptotically flat P-L or NS solution there are no event horizons, 
thus A(p) has a constant sign for c < p < oo and because A(p) > at spatial infinity so it is for c < p < oo. Only 
for P-L solutions, the functions A(p), V{p) & <p{p) are finite everywhere for c < p < oo. 

Some static, spherically symmetric, asymptotically flat solutions to Eqs Q to (H) have been derived in [2, 
4, 6], some of which represent different configurations of singular BHs. In Ref [6], Bronnikov and Shikin 
constructed BH, NS and P-L solutions. For the latter solution, their function A (Eq (B.3) of Ref [6]) does not 
approach 1 at spatial infinity (p ^ oo); rather, it approaches the limit 1 +2(tanhc— l)sech^c, where c is their 
parametric constant. Thus, their "P-L" solution is asymptotically flat only in the limit c — > oo, corresponding to 
the trivial case r = p with zero mass and constant functions A = 1, (j) = \/37r/\/8 & V = 0: this is Minkowski 
space time. Recently, Bronnikov and Chernakova [2] constructed a two-parameter asymptotically flat BH char- 
acterized by its negative potential V and where V & <p diverge near the center and converge as p — > oo, while A 
is finite everywhere. Very recently, Nikonov et al [4] have obtained a two-parameter BHs and NSs as well as 
a one-parameter P-L solutions. For the latter, all the three functions A > 0, & V are finite everywhere, V has 
both signs and is an incomplete elliptic integral of the first kind. 

Most of the solutions constructed previously [2-4, 6] were derived using the inverse problem method. The 
method consists in choosing a suitable function r(p) then solving Eqs Q to dSJl for A, V & (Eq ^ is auto- 
matically satisfied). There are physically many interesting exact solutions to Eqs (l5]l to ([H) depending on the 
particular form of r(p). The behavior of the solutions at spatial infinity is somehow defined by the boundary 
condition requiring that the solutions be asymptotically flat 

A(p) = l-2m/p + C>(l/p^),r~p(p^oo) and lim / = 1 , lim r" = 0, (9) 

where m is the Schwarzschild mass. At the center there is no specific boundary condition, this is why there is 
a variety of possible asymptotically flat solutions, however, grouped into three categories: different (singular) 
BH solutions, NS and P-L solutions. To be more specific, the BH set (A having both signs) include solutions 
where V is finite everywhere (this paper) as well as solutions with divergent V [2] ; the same applies to the NS 
set where A > (this paper). The solutions discussed in [2,4,6] and those discussed in this paper have different 
boundary conditions at the center. 

It is an interesting issue to investigate the behavior of the generic solutions near the center to determine the 
criteria by which the solutions emerge. The solutions discussed in [2, 4, 6] have been derived by choosing a 
special f unctioiH r(p), however, they lacked a systematic treatment and many questions remained unanswered. 
For instance, we do not know a) whether it is possible to obtain a BH (resp. NS) where both V & are finite 

^The method developed in [4] for obtaining the solutions to Eqs l|5]( to ([Sjl looks somehow different from that developed in [2,6], 
however, it is equivalent. In fact, in [4] the function r(p) has been chosen such that dr/dp = (r"* +4)/(r'* +4 — 4a'*), where a is their 
parametric constant. 
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everywhere, b) what conditions are required for obtaining, say, a BH or a P-L solution, c) under which conditions 
the potential V has both signs, etc. 

In section m we derive new integral formulas for A, V & <p and proof five corollaries and a theorem on 
selection criteria for obtaining predefined static, spherically symmetric, asymptotically flat solutions to Eqs ([5]l 
to (H]). First we deal with the most generic case then focus on a constrained case with minimum conditions. In 
section [3l the behaviors of A, V & ^ near the center and at spatial infinity are discussed, instances of solutions 
of particular interest are provided and two energy conditions are analyzed leading to two further corollaries. 
In section |4l we derive further solutions by the inverse problem method. Two appendices have been added for 
completeness. 



2 Criteria for obtaining predefined asymptotically flat solutions 

Eq ^ is readily brought to the form [A(r^)' — r^A']' = 2 which is integrated by 

r'^A'-{r^yA = -2p+K, (10) 
where K is constant of integration. Integrating (ITOl) leads to 

Aip)=Cr^{p)+^f^dzyHp), 

where C is a second constant of integration. With the boundary conditions we must have C = and K = 6m. 
Hence, 

Using [A{r'^y - r^A']' = 2 into Q, we bring V to the form 

, , (r^A'Y 1 (Ar/)' 
np) = -V = ;,-V. (12) 

Substituting ([TTI ) into the one or the other expression of V, one obtains 

V{p) = -,- ^ - (3/^ + rr") . (13) 



Finally, without loss of generality we take the negative sign in ^ = — \J —2r" jr) and introduce two arbitrary 
constants po > c & 0(po), we obtain 

(pip) - (piPo) = - r V-^r"{z)/r{z)dz. (14) 

"'Po 

Throughout the remaining part of the paper, when the limits of A, V & <p at the center c exist, we define 
A(c), V{c) & (j){c) by limp^c+^(p) ^^p^c+^iP) & linip^c+ 0(P)' respectively. This definition extends to 
their derivatives at the center. 



2.1 Generic case 

Notice that Eqs ([TT]) . ([T3] ) & (fT4l ) have been derived without referring to any extra condition (except the boundary 
conditions (O and, in order to integrate (ITOl) . we have assumed implicitly that /(p) exists (r(p) differentiable)), 
thus the following first two corollaries are generic. 

Consider the case where c > 3m. Since 3m < c < p < oo we have 3m — p < 0, thus the integral in (fTTI) is 
positive leading to A > 0. Hence, 
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Corollary 1 If the mass of the spherically symmetric, asymptotically flat solution to Eqs ^ to ^ is such that 
c > 3/M where c is the scaling constant, then the configuration is either a NS or a P-L solution. A negative-masse, 
spherically symmetric, asymptotically flat solution is either a NS or a P-L solution. 

Now, let us look for the conditions that make (the limit of) A finite at c. Writing = + /^f in ([TT]) . 
where u is any fixed number in the open domain (c, oo). The improper integral dz converges by the limit 

comparison test. In fact, let /(p) = \3m — p\/r^ & h{p) = 1/p^, since Yimp^^f/h = 1 and since J^h{z)Az 
converges for c < m < oo, then the integral of / is convergent. With that said, the value of A at c is the following 
limit 

limA p = 2 lim , !f . ^ - , (15) 
' lim^. (16) 



/(c) p^c+ r(p) 

The conditions /' < (Eq ©land hmp^oo r' = 1 (Eq I©) imply 1 < r'(p) < oo for c < p < oo. Now, with r{c) = 
0, the limit in ([T6l ) exists onl>|^if limp^f.+ (3m — p) = 3m — c = 0. In this case A (c) = limp^c+ A(p) = l//^(c). 
If c > 3m, limp^f+ A(p) = +oo and the solution is a NS by corollary 1. If c < 3m, limp^£.+ ^(p) = ~°° the 
solution is a BH since A (p) has both signs. In this case (c < 3m), A(p) is manifestly positive for 3m < p < oo. 
For c< p < 3m, we can write the integral in ([TT]) as j£ = + f-^^^^ where f^ ■■■ = constant > and J^^ • • • < 0. 
As p decreases from 3m to c, the second integral J^^^^ ■ ■ ■ decreases monotonically from to — oo (since we found 
limp^c+^(p) = — oo). Consequently, the sum of the two integrals and A(p) vanish at one and only one point ph 
(event horizon) between c and 3m: c < ph < 3m. Thus for c < 3m, ph is a simple root of A(p) =0 (all horizons 
are simple). 

Corollary 2 If r{p) is differentiable at c, a necessary condition to find a spherically symmetric, asymptotically 
flat P-L solution to Eqs ^ to ^ is c = 3m and in this case A(c) = l/r'^(c). In the other cases, a) c > 3m 
including negative masses, the solution is a spherically symmetric, asymptotically flat NS with limp^c+ A(p) = 
-foo or b) c < 3m, the solution is a spherically symmetric, asymptotically flat BH with limp^c+'^(p) = ~°° ^f^d 
the event horizon c < p^ < 3m. 

The condition c = 3m in the above statement for obtaining P-L solutions being necessary, one may derive 
solutions with c = 3m which are not P-L solutions. By corollary 1, these solutions must be NSs. Thus, corollary 
2 does not rule out NSs with c = 3m. Solutions with c = 3m and limp^c+ ^(p) divergent are indeed NSs. 

We shall see in the next subsection that, for the conditions on which corollary 2 rests, A(p) is finite for 
< p < oo. Since the condition c = 3m ensures finiteness of A at c, thus by corollary 2 we have solved the 
eigenvalue problem consisting in determining an everywhere finite solution A(p) to Eq ^ with the "eigenvalue" 
m = c/3. 

Repeating the steps ([T5] ) & ([T6l ) for limp^(.+ (rA), we obtain 

lim [r(p)A(p)]= 2 lim /ff = " P) = ^^?77T^ ' ^7) 

p^c+ p^c+ l/r^(p) 3r[c) p^c+ 3r[c) 

Hence for c 7^ 3m, A behaves near the center as 



■'An instance where A(c) is finite with {3m — c) arbitrary is r = \/ — (?-■ However, for this function r' & r" are divergent at i 
consequently V too diverges at c, thus there is no P-L solution. 
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2.2 Constrained case with minimum conditions 

In addition to the boundary conditions Q and differentiability condition of r{p), we impose further (minimum 
set of) conditions on r(p) and its first derivatives in order to derive further selection criteria. We will do that 
without losing much generality. Thus, throughout the rest of this section, we assume that r(p) is four times 
differentiable in the domain c < p < oo, which means that there r(p) and the derivatives r'(p) to r^^\p) are 
smooth functions (continuous) and finite 

r(p), r'(p) — > r^'^\p) : exist and continuous onc<p<oo (19) 
< r(p) < oo, 1< r'(p) < oo, -oo < r"(p) < 0. (20) 

Notice that /(p) > 1; particularly, we have 1 < /(c) < oo. In Part 1 of Appendix A, however, we will see 
that the third and fourth derivatives of r(p) need not exist on the whole domain p € [c,oo); rather, on a small 
half-open interval containing c. 

The conditions (fT9l ) & (l20l ) do not ensure existence of Taylor series expansions or others for r(p ) at the 
center c since the fifth and higher order derivatives of r(p) are not assumed to exist there and elsewhere. The 
same is true for A(p), V{p) & <p{p), which are related to r(p) via Eqs (fTTI) . (fT3] ) & ([14)) . For this reason the 
method we are considering in this subsection (as was the case in the previous one) and Appendix A is more 
general, different from and even advantageous than series expansion methods. 

We start by proving the finiteness ofV{p) & 0(p) at spatial infinity (that of A(p) is imposed in Q). Us- 
ing ^ in ( fT3] ) & (fT4l ) one obtains limp^co V = & limp^oo = constant. The result for V is obvious and to 
derive that for (p we take the limit as p ^ oo of both sides of ([T4l) 



lim(^»(p) = 0(po)-lim r V-2r"(z)/r(z)dz. (21) 



We can always bring r(p) to the form r(p) = p — s{p) such that Hmp^oo^' = 0, limp^^o^' = & limp^oo^" = 0. 
More precisely, as p ^ oo, ^ < a^/p" where n is some positive numbei^. Hence, as p ^ oo, —2r" (z) / r{z) ~ 
a^j2n{n + l)/2("+3)/2 and the integral in the r.h.s of (|2T]) converges to some constant. 

For c<p<oo(0<r< oo), the improper integral in (fTTI) converges by the limit comparison test, which 
means that A is finite and continuous for c < p < oo. A{p), r'{p) & r"{p) being finite functions in the domain 
c < p < oo, Eq ([T3] ) implies that V remains finite for c < p < oo. Finally, since the function 2r"(p)/r(p) is 
continuous on any closed interval included in the open domain (c,oo), the integral in ([T4l) exists for po > c, which 
means that is finite and continuous for c < p < oo. We have thus shown that the three functions A , V & are 
finite and continuous for c < p < oo; the unique point where the functions may diverge is the center c. 

Notice that the continuity and regularity of the functions A, V & for c < p < oo is independent of the 
nature of the solution (BH, NS or P-L solution). Thus 

Corollary 3 Ifr{p) is constrained by (1191 ) & (1201 ). no singularity at finite r (resp. p) occurs in the open domain 
(0,oo) (resp. (c,oo)j whatever the nature of the solution is. 

Now, we look for the conditions that make the (limits of the) functions V & finite at c (that of A has been 
established in corollary 2). The behavior of V near the center is derived from ([T3] ) as follows. 

\mi[r^{p)V{p)] = 1 - lim [(rA)r"]+ lim [rr(p)] (22) 

3r(c) p^c+ 



*In the worst case, s{p) behaves as In '^p/p ()f > 0). Thus as p ^ oo, i- < 1 / ^/p {n=\ /2). 
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where we have used ([17] ) to evaluate the first limit in (l22l) . In (|23] ). T{p) is defined by 

T{p) = 2(p-3m)-3ArV ^ 
r 

which is by ([TT] ) an indeterminate form 0/0 as p ^ c+. Using Eq ([TOl l with = 6/m, it is straightforward to show 
that T = -{A'r+Ar'). Thus 

,. 2(p-3m)-3Ar'r 
lim r = lim — 



9 "^Arr" 

3 lim (A V+ A/)- lim 



r'{c) p^c+ p^c+ r' 

2 , ^ 2(c-3m)/'(c) 

■ + 3 hm r - , ^ ^ . (24) 



/(c) p^c+ r'2(c) 

Solving (l24l ) for limp^c+ 7" then substituting back into (1231 ) yields limp^^.+ (r2y) = (c — 3m)r"(c)/(3r'(c)). 
Thus, for c / 3m we have 

(c-3m)/'(c) 1 , 

^(p^" 3/(c) 

Corollary 4 If r{p) is two times dijferentiable with \ <r'{c) <o° (r{p) need not satisfy all the conditions (I19l l 
& (1201 )) a«<i ifr"{c) ^ 0, f/j^n as p ^ a) the function V (resp. A) for a NS solution with c > 3m (corollary 
2) approaches — oo as —l/r^ (resp. +oo as \/r) and b) the function V (resp. A) for a BH solution (c < 3m by 
corollary 2) approaches +oo as l/r^ (resp. — oo as — l/r). 

Concerning extreme NSs with c = 3m, we have the following result. 

Corollary 5 If r{p) is two times dijferentiable with 1 < r'{c) < oo (r{p) need not satisfy all the conditions (1191 ) 
& (l20l )) and if r"{c) / 0, then as p ^ c^ the function V (resp. A) for the NS with c = 3m approaches — oo as 
— l/r (resp. converges to \/r''^{c) by corollary 2). 



Now, we will look for conditions to have a P-L solution and provide a proof for corollary 5 (Eqs (1271) . 
& (l30l)). By corollary 2, P-L solutions exist only if c = 3m. 

Let c = 3m. In Part 1 of Appendix A, we will prove the general result that if c = 3m then limp^c+'^' 
converges. With K = 6m = 2c, Eq ( flOl ) reads A'(p) = [2r'rA + 2{c — p)]/r'^. Taking the hmit of the previous 
formula as p ^ c"*", replacing limp^f+ A' by A'(c) and using limp^^+A(p) =A(c) = l/r'2(c) we obtain 

2r'rA + 2(c-p) 
A c = hm 2^ - 

, /VA + r'rA'+(/2A-l) 
A (c) = hm ^ 

/. ^ ''"(c) „ ^ 1 (/^A-l) 
r^(c) r(c) p^c+ r 

Dropping A'(c) from both sides of (1251 ). the remaining equation reads 

r"{c) (/2a -1) 

r^(c) p^c+ r 

/'(c) 2//'A + /2a' 2/'(c) , . 

/2(c) p^c+ / /2(c) p^c+ 
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Solving the last equation for A'(c) we obtain 

A'(c)=^lim^A'(p) = -^. (27) 

Eq ([D reads 

rA" = 2Ar" + 2{r'^A-l)/r. (28) 
Using (ED and limp^,,+ A(p) = A(c) = l/r'^{c) in (EH) we obtain 

2r"(c) r"(c) 

lim(rA") = ^-2-^=0. (29) 

The first equation in ([T2l) reads 

V,p, = -?:^. ,30) 

2r 

Notice that by ([271) & (l29ll, limp^,+ {Ir'A' + rA") = -6r"(c) //^(c), thus if r"(c) / (r"(c) < 0) then limp^c+ 
V (p ) approaches — oo as — 1 /r; this establishes a proof for corollary 5 . In this case V &.A behave near the center 
as 

Now, we further assume that r"{c) = (with c = 3m), the limit of both sides of (l30l ) as p — > c+ leads to 

A' A" 

lim V = -r{c) lim lim — (31) 

p^c+ p^(:+ r p^c+ 2 

= -/(c) lim ^-Y,ra^ = -l lim A" . (32) 

Since r"(c) = this implies A'{c) = by (|27] ) making thus the first limit in the r.h.s of (|3T]) an indeterminate 
form 0/0. Next, consider the second term in (fT2l ) 

l-AV-Ar'2-Arr" 
lim V = lim r > - 

p^c+ p^c+ 

/ ,A' A'r" 3A /' A/" A"\ 
p^c+ V?" r 2 r 2r 2 J 

where the second limit in (l33l) vanishes by the imposed condition: r"{c) = 0. The first and third limits are 
indeterminate forms 0/0, leading to upon further differentiating 

3 2r"'ic^ 
lim V = -- hm A" - „;/ . (34) 

p^c+ 2 p^c+ r'^{c) 

Eqs (l32l) & (iMl) result in 

lim V = - ^ hm A" & r'"(c) =0. (35) 

In Part 1 of Appendix A we will prove the general result that if c = 3m & \5r"^{c) — Ar' {c)r"'{c) = then 
limp^c+^" (=^"(c)) converges and provide a formula for A"(c). In our case, since r"{c) = & r"'{c) = 0, the 
latter condition is satisfied and the desired formula for A"{c) will have the form (see Eq (IA.9I )) 

A"(c) =2 lim {r^A-\)/r^ = L, ifc = 3m, r"(c) = 0, &r"(c) =0. (36) 

p^c+ 
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Table 1: Static, spherically symmetric, asymptotically flat solutions to Eqs Q to (O. Column 1 from the left shows the 
statements concerning the behavior of the derivatives of r(p ) and their values at the center p = c (r = 0). At spatial infinity 
p =°°, the boundary conditions (|9]l are assumed for all solutions. In Columns 2, 3 & 4 the behaviors of A(p), V{p)&(j){p), 
when available, are given in the limit p ^ c+ (r ^ 0+) and only the dominant terms are shown. Notation & nomenclature: 
A0 = ^(p) — 0(c), x = p—c, "NC" for "Necessary Condition", "SC" for "Sufficient Condition". In the expressions of V 
depending on Inx, we assume existence of A^'^ (p ) as stated in Eqs ( |38] | to i 



CONDITION(S): 

At p = oo the conditions 
^ are assumed 



Particle-Like Solution 
A{p)>Q 



Naked Singularity 



A(p)>0 



Black Hole 
A(p) has both signs 



Corollary 2: 
r'(p) exists 



c = 3m (NC): 



a) m < c/3 (SC): 



3r'(c) rip) 
b) c = 3m: Possible 



m > c/3 (SC): 



3r'(c) rip) 



Corollaries 4 & 5: 
r'(p) & r"(p) exist, 

r"(c)^0 



a) m < c/3 (SC): 



Vip) 



3r'(c) r(p) 
c-3m)r"(c) 1 



m > c/3 (SC): 
^(p).,2(c-3m) 1 



b) c = 3m: 

A(p): 



3r'(c) r2(p) 



V(P) 



3r'(c) r(p) 
(c-3m)r"(c) 1 

37(c) r2(p) 



1 



y(p) 



r'2(c) 
3r"(c) 1 



r'2| 



Theorem: 
r'(p),r"(p),r"'(p) 

& rW(p) exist, 
r"(c)=0, r"'(c)=0 



c = 3m (SC): 

Mp)^:j^&a'{c) = o 



Vip): 



r'2(c) 
-3L 



m < c/3 (SC): 
^(p).,2(c-3m) 1 



3r'(c) r(p) 



m > c/3 (SC): 
^(p).,2(c-3m) 1 



, , (3m-c)rW(c) 
y(p) ~ -^^ . J . ^ Injc 



3r'(c) rip) 



A0 ~ - 



2 /-r(4)(c) 



r'ic) 



.3/2 



A0 



r'3(c) 

-rW(c) 



y(p)^(^-.f)-!"(-)in. 



,3/2 



A0 



r(4) 



(C) 



r'(c) 



,3/2 



Notice that 
0. 

Eq ([HI) implies 0'(c 
(-2r'7r)]i/2 = -[lim 



is derived from (1281 ) after dropping the first term, which is zero: limp^f.+ r" /r = limp^c+ /" // ■ 



limr 



limn 



,+ y=2?Vr=-[lim^ 



ip^c+ ("2/"//)]^/^ = (by (l35l)). thus is differentiable at c and consequently it is 
continuous there, which means that 0(c) exists and is finite. Finally, we managed to make the three function 
A, V &<p finite at c, thus finite everywhere in the domain c < p <°° under the following conditions. 



Theorem 1 Suppose r{p) is four times differentiable at every point of the domain c < p < oo and satisfies the 
conditions ( [T9l ) & ( |20l ) (/" & need not exist on the whole domain p £ [c,°°); rather, on a small half-open 
interval containing c). Ifr"{c) = 0, r"'{c) = Q & c = 3m, any spherically symmetric, asymptotically fiat solution 
to Eqs ^ to ^ is a one parameter P-L configuration withA{c) = l/r'^(c), A'(c) = & V{c) = — 3A"(c)/2. 



Table \\\ summarizes the results of the above discussions and includes some results from the following sec- 



tion. 



Most of the solutions constructed previously [2-4,6] were derived using the inverse problem method, which 
we are aiming to use in section [4] The method consists in choosing a suitable function r{p) then using Eqs ([TT]) . 
([n) & CD) to obtain A, V 
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From the above discussions (five corollaries and a theorem), we have learnt how to set conditions on r(p) 
and to fix the sign of (c — 3m) in order to obtain a predefined final solution. For instance, for r(p ) given by 

r(p) = ^{p-c)--^ ^^^^ ,p>c>0,n (integer) > 2, (37) 
n—l n — I p" ^ 

(where r"{p) = -nc^ {p - c)"-^ / p"+^ < 0) if we choose n > 4, then r(p ) satisfies all the requirements of the 
theorem but one; if we further choose c = 3m in (fTTI ). then Eqs ([TT]) . ([T3l) & ([T4l ) lead to a one-parameter P-L 
solution for each value of « > 4. If, instead, we choose c > 3m or c < 3m the solution is by corollary 2 a 
two-parameter NS or a BH, respectively, for each value of « > 4. 

Notice how crucial is the condition r"'{c) = for obtaining P-L solutions for the set of the continuous 
everywhere four times differentiable functions r(p) satisfying the conditions ( [T9l ) & (l20l) . For n = 3, this 
condition fails (r"(c) = & r"'{c) = —3/c'^) and the above formula ( [37] ) gives rise to no P-L solution, whatever 
the value of (3m — c) is; even when c = 3m, the potential V still diverges at c as r"'(c)ln(p — c) (see Eq (|4T])). 
Similar conclusion holds for « = 2. 

In the following section, we will derive further properties of the solutions at both the center and spatial 
infinity for more restrictive conditions on r(p). We will also discuss the null and weak energy conditions. 

3 Generic properties and energy conditions 

Here we study the properties of the asymptotically flat solutions (in their generic forms) by investigating their 
behavior near the center and at spatial infinity. We also examine the weak energy condition (WEC) at these 
end-points and the null energy condition (NEC) on the whole range of p. 

3.1 The center vs. spatial infinity 

In this section, we rather assume that r(p) is five times differentiable: 

r(p), r'(p) —>■ r^^^(p) : exist and continuous on c < p < oo (38) 
< r(p) < oo, l< r'(p) < oo, -oo < r"(p) < 0. (39) 

The behavior of A near the center is discussed in 1) Part 1 of Appendix A for r(p) satisfying (fT9l) & (l20l ) 
with c = 3m and 2) Part 2 of Appendix A for r(p ) satisfying (l38l) & (l39l) with c / 3m. 

If r(p) satisfies (1381 ) & (|39l). A. V & <p take the following forms near the center p = c (r = 0), where A"{c) 
is determined in Part 2 of Appendix A by ( IA.17I ) 

2(c-3m) 1 1 4{c-3m)r"{c) F 
3/2 (c) x^THc) 3?Hc) ^V^c)^ 

+^xHnx + ^^ + ..., (40) 
r ^ 2 



V 



(c-3m)r"(c) 1 \Q{c-3m)r"^[c)-r'{c)[9r"{c)+A{c-3my"{c)] 1 
72 



3/3 (c) 



3/4(0) 



3A"(c) \A{c-3m)r"\c) 10/" (c) 157(c - 3m)/'(c)/"(c 



2 /5(c) 
3(c-3m)rW(c) 14/'^' 



4/3 (c) 



■ + ■ 



3/3 (c) 
3G 



+ ■ 



/4(c) /5(c) 



lnx + 



H 



18/4(c) 



-xlnx + ■ 



(41) 
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where x = p — c and F, G & H are polynomials of {r'{c),r"{c),r"'{c),A'^^ , (c — 3m)). The polynomials F & G 
are given by 

57 (c - 3m)r"2 (c) _ r'{c) [54/' (c) + 20(c - 3m)r"'{c)] 
F = 



18 

G = " — ^-^-'^''^ +5(c-3m)/'^(c) 



4r'"(c)/2(c)-15r"2(c)r'(c) , 



10(c-3m)/'(c)/''(c)/(c) (c-3m)rW(c)/^(c) 
3 ^ 3 ■ 

From now on and throughout the rest of this paper, we choose po = °° and (j) {o°) = (precisely limp^oo (p ) — ) 
0+) in ([T4l ). The function r(p) being concave down on its range, the first nonzero derivative r^"\c), 2 <n <5, 
is negative (r'(c) > 1 always). The behavior of 0(p) depends of the first non-zero A"\c): 



Hp) = Hc)-VsJ^^V^+---,ifr"{c)^o, 



Hp) = Hc)-V2J^^x + ...,ifr"{c)=0,r"'{c)^0, (42) 



Hp) = 0(c)-^y'^J^xV2 + ...,if/'(c)=O,/''(c)=O,rW(c)/O, 

where 0(c) > (H{p) < and we have chosen limp^oo0(p) 0+). 

The behavior at spatial infinity is somehow defined by the boundary condition This condition Q does 
not completely fix the behavior of A, V & (j) at spatial infinity. We fix that by requiring r(p) of the form 
r (p ) = p + ao + Lr= li'^n/p)- The constant ao can be absorbed in the first term upon performing the translatioiil 
p ^ p + fl;o- Thus, at spatial infinity we assume r(p) of the form 

CO 

Kp)=P + LK/p)- (43) 

This is the simplest form to which one may bring r(p) to. The solutions discussed in [2,4,6] have the same 
expansion as in ( |43] ). 

Upon substituting (l43l) into ([TT]) then into ([T3l) . one obtains 



, , 2m 2{2ma\+a2) -\-2a-i , , 



5p-' p 

Notice that, since the function r(p) is concave down on its range, the first nonzero constant a( in (1431 ) is 



^Such a translation modifies the value of the new p at r = 0, which is now c + ao, without modifying the scaling constant (= c). In 
Appendix A, we perform a similar translation and work with a radial coordinate x whose value at r = is instead of c. Such a translation 
is no trouble, only some equations, like Jilt , change the form (compare with l lA. 12l (). however, the derivations and conclusions are the 
same. 
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negative. Accordingly, the series expansions for 0, which are derived from (fT4l l or read 



Hp) 



1^ —a\ Saj 



2+0(l/p^), ifai /O 



p 2^/^p2 



+ 0(l/p^/2) 



if ai = 0, a2 / 



^p3/2 57^pV2 



0(P) 



^/—6a3 5V2a4 



- + O(l/p4),ifai=a2 = 0,fl3/0. 



p2 3V3y^p 



One draws the following conclusions concerning the behavior of V, ^ & A as p ^ oo. Although V ^ as 
p ^ oo, it may approach zero from above or from below by adjusting the parameters £?„. For instance if ai ^ 0, 
02 & aim > 3(32/4, the potential is positive at spatial infinity, which means the weak energy condition 
(WEC) is satisfied. Since in this case ai < 0, the previous condition on m is satisfied if m < 3a2/(4fl'i): the 
mass of the solution is bounded from above by 3a2/(4ai). Particular solutions where the mass is bounded from 
above have been obtained in [4]. Conversely, if the mass is bounded from blow m > 3a2/(4ai ), the potential is 
negative at spatial infinity. We will show in the following subsection that in this case the WEC is also satisfied. 
If, otherwise ai = a2 = but a^, / 0, the potential is always positive at spatial infinity (V m) since in this case 
as < an again the WEC is satisfied. The scalar ^, which also vanishes at spatial infinity, has always the same 
sign as p — > oo. Finally, for m>0,A always approaches its limit 1 from below. 

Corollary 6 At spatial infinity, we assume that r{p) — p admits an expansion in powers of l/p given by (1431 ). 
Then ifa\ (ai <Q) & a2^0, there is a critical mass mc = 3a2/(4ai ) beyond which (m> m^) the potential 
is negative at spatial infinity; if, instead, m < m^ the potential is positive. For m = m^. the potential may have 
any sign. In all cases, the WEC is satisfied at spatial infinity. 

In the paragraph preceding corollary 3 we have mentioned that the finiteness of A, V & for c < p < oo is 
independent of the nature of the solution. A first important and yet general similar result from the set of Eqs (l40l) . 
dHJ & (|42)) is that is always finite at the center, provided conditions (l38l) & (l39l) are satisfied. Although A(c) 
is finite for c = 3m, as seen from Eq ( |40l ). it is obvious from Eq (|4TI ) that the condition c = 3m alone does not 
ensure existence of P-L solutions for, in this case, V diverges near the center as 3r"{c)/[r'^ {c)x\ — > — oo and the 
solution is a NS. Similarly, the conditions r"{c) = & c = 3m do not ensure existence of P-L solutions for V 
diverges near the center as —4r"'{c)lnx/r'^{c) —>■ —°° and the solution is again a NS. 

Other important conclusions derived from Eqs ( |40l ) & (|4T]) . not discussed yet in the literature, follow. Under 
constraints (EH) & (Ull: 

a) The field Eqs (l5]l to ([8]l do not admit BH solutions (c < 3m) where the metric A is finite everywhere. However, 
one can derive a BH solution where both V & (j) are finite everywhere by choosing r"{c) = r"' {c) = r'^'^\c) = 0; 
in this case, A diverges at the center as 2(c — 3m) / [3r'^ {c)x\ —°° ( (l40l ) or ([TS])). One may call such a solution 
a "semi-regular" BH since all matter fields are finite. We have a set of such two-parametric solutions given 
by (|37] ) provided we take « > 5 to ensure that r"{c) = /"{c) = rW(c) = and in (dB we choose c < 3m. The 
simplest of these solutions is 



In Ref [2], another "semi -regular" black hole solution with A finite everywhere has been constructed, how- 
ever, by choosing a function r(p) not differentiable at the center for which limp^c+ f'{p) =°° (in our solutions 



/X 5 1 (P-' 

Kp) = 4(P-c)-4^ 



p >c>0 



c < 3m . 



(46) 
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1 < r'{c) < oo). This choice of r(p) caused V & to diverge at the center. 

b) Similarly, if c > 3m the field Eqs Q to dSJl do not admit NS solutions where A is finite everywhere. However, 
if c > 3m one can obtain a NS where both V & are finite everywhere by choosing r"{c) = r"'{c) = r^^^ (c) = 0; 
in this case, A diverges at the center as 2(c — 3m)/[3/^(c)A;] — > +oo ( (|40l ) or ([TS])). We have a set of such two- 
parametric solutions given by (|37] ) provided we take n > 5 to ensure that r"{c) = r"'{c) = (c) = and in (HB 
we choose c > 3m. The simplest of these solutions is 

r(p) = ^(p-c)-^^^^,p>c>3m>0. (47) 

c) For any NS with c = 3m, A is finite everywhere and V diverges at the center as —4r"'{c) Inx/r'^ (c) —>■ —oo if 
(/'(c) = 0&r'"(c) / 0)oras3/'(c)/[r'3(c);c] ^-ooif r"(c) /O. Asolutionof thecase (r"(c) = & r'"(c) / 0) 
is (Eq ^ with n = 3) 

r{p) = l{p-c)-^^^^^, p>c = 3m>0, (48) 
and a solution of the case r"{c) / is (Eq (|37] ) with « = 2) is 

(p-c)^ 

r{p)=2{p-c)-— p>c = 3m>0. (49) 

P 

3.2 Null and weak energy conditions 

For the case of self-gravitating real scalar field, the NEC and WEC are the mostly discussed conditions in the 
literature [5,9]. If & are any future-directed null {k^k^ = 0) and timeUke {u^u^ = 1) four- vectors, the 
NEC and WEC state 

NEC: r'^v^^,^" >0, WEC: r^^vM^"" >0- (50) 

It was established in [5] that the conditions which rule out wormholes, horns and flux tubes rest on the 
validity of the NEC. In this section we will show that (P-L, NS or BH) solutions to Eqs ([Sll to ([8]) do not violate 
the NEC. Hence, the consideration of the NEC does not allow for further selection criteria. 

Let £ denotes the sign of A such that £^ = 1 & £A = |A| (£ = — 1 only inside a BH). A tetrad frame associated 
with the metric ^ is given by 

eo^ = (£x/H, 0,0,0), ei^ = (0, -£/v^,0,0) , = (0,0, -r,0), 

e3M = (0,0,0, -rsine), ^'o'^ = (l/y^, 0,0,0) , ej'' = (0, /H,0,0) , (51) 
= (0,0, l/r,0) , e^^ = (0,0,0, l/(rsin d)) . 

Future-directed null {k^) and four- velocity {u^) vectors take the forms 

ki" = eo^ + siei^ + S2e2^ + s^e^^ , with e{l - s{) - s^- sl = , (52) 
=N{eo^ +Siei^ +S2e2^ +Siei^), with = l/yfiO^^^^yr^TT^. (53) 

For the static, spherically symmetric solutions to Eqs ^ to ([S]), the fields A, V & depend only on p. 
Eqs © & ® lead then to T', = T^e = T'P<p = V + {A^''^/l), - {A^'^/2) and r'^v = if / V. 

Hence 

T%k^k'' = esj{T\-TPp)=si\A\(j)'\ 
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Since > (the r.h.s is zero if = (with £ = 1) or/and r"{c) = 0), solutions to Eqs ([5]) to (Ull do not 

violate the NEC. Notice that the NEC is fulfilled even if the field Eqs © to ([HI) are not satisfied. 
Similarly, we obtain 

T^yu^.u" =N\e -si- sl)T't - es]TPp] 

=N^{[e{\ - sj) -si- sj]V + [l+sj- esl - £^i]l^l0'V2} • (54) 



The number in the first square parenthesis of (l54l) is positive by (1531 ) (A'^ real). This leads to 1 +>yj — esl — es^ > 
2s J if £ = 1 or 1 + — £^'2 — £5'3 > + 1 if £ = — 1 . Thus, at the points where V >0, the r.h.s of (l54l) is positive 
and the WEC is satisfied. 

The sign of the r.h.s of (l54l ) cannot be determined without knowing the solution to Eqs ^ to ([8]). Since the 
functions A, V & (p' are related to r(p) and its derivatives via Eqs ([TT]) . Q & Q, it suffices to know the shape 
of r(p) to completely determine the sign of the r.h.s of (l54l ). In this paper we have chosen r(p) by imposing the 
conditions Q & (1431 ) at spatial infinity and the conditions (1381 ) & ( [39l ) at the center. 

Consider r(p) of the form (1431 ) and let p*^ & pp be the mass density and radial pressure due to the field (p. 
Expressing A, V & (j)' in powers of 1/p, we obtain 



'(£ + 3)p^+3 ^ P {£ + 3)pM 



(at spatial infinity e = I & I — sj — S2 — si^ > 0). Here is the first non-zero constant a„ in (1431) : ai < 0. Thus, 
the dominant term in the r.h.s of ( [55l ) is positive. 

Corollary 7 a) The spherically symmetric, asymptotically flat solutions to Eqs ^to^do not violate the NEC. 
b) At the points where V >0, the WEC is satisfied and the mass density is positive, c) If the fimction r{p) — p 
admits an expansion in powers of 1 /p, then at spatial infinity the spherically symmetric, asymptotically flat 
solutions to Eqs ^ to ([8]l, whether V > Q orV <Q, do not violate the WEC and the mass density and radial 
pressure are positive. 

The consideration of the WEC does not allow for further selection criteria at least at spatial infinity; this 
is not the case near the center where it is generally possible to select within the same family (P-L, NS or BH 
solution) those solutions which do not violate the WEC, as will be established shortly below. For that purpose, 
we focus on the case of a BH solution then tabulate the results for the other solutions. 

We assume the conditions (l38]l & dlHl hold and use Eqs (gOll & (|4T]|. For a BH solution c - 3ni< 0, £ < & 
s\>l+sl + s]. If r"(c) / (/'(c) < 0), in the hmit x ^ 0+ Eq (m leads to T^yu^u" oc (c - 3m)r"(c) /x^, so 
the WEC is satisfied. If r"{c) = 0, T'^yU^u" oc (c - 3m)r"'{c) /x so that the WEC is satisfied if /"(c) < 0. Now, 
if /'(c) = & /"(c) = 0, we have T^yu^u" oc -(c - 3m)r('^)(c)ln;c and the WEC is satisfied if A'^^c) < 0. 
Finally, if /'(c) = 0, /"(c) = & r'^) (c) = 0, T^yU^u" oc -L and the WEC is satisfied if L < 0. 

Table |2] shows the conditions for which the WEC is satisfied near the center. There is no case of validity 
of the WEC for a NS with c - 3m > & /'(c) = (case lb) in Table IS): the WEC is violated V /"(c) 7^ 0. In 
fact, under these conditions we have T^yU^u"^ oc (c - 3m) (1 - 3s? -si- sl)r"'{c) /x. The WEC is violated if 
X-'isl-sl-s'l & r"'{c) have opposite signs. Since for a NS £ = 1, < Sj < 1 & < sf + ^2 + ■^'s < 1, this 
implies —2s\ < 1 — ?)s\ — si — s\ < \ — ls\ < 1, thus 1 — 3s\ -sl-s^ may have both signs. For instance, for a 
four- velocity vector or an observer with si = 1/2 & ^2 = ■^'3 = 0, 1 — 3sj — ^2 — ■^•3 = 1/4 and the WEC is violated 
if /"(c) < 0; for an other observer with si = 1 /\/2 & S2 = S3 = 0, \ -3sj - si - sj = -1/2 and the WEC is 
violated if /"(c) > 0. Similarly, there is no case of validity of the WEC for a NS with c — 3m = & /'(c) < 
(case 2a) in Table H]). 



14 



Table 2: Static, spherically symmetric, asymptotically flat solutions to Eqs (|5]l to ((SJl. Validity of the WEC near the center 
for P-L, NS & BH solutions. r(p) is subject to conditions (l38) & ( l39] l. 



Parttci F-T tkf 


Naked Singularity: 


Rt ack Hot F' 




Solution: c - 3ot 0, 


c — 3m > 




c — 3m < 




r"(c)==0, r"'(c)=0 












1) c-3m > 0: 










a)r"(c) <0 




a) r"(c) < 










b)r"(c)=0&r"'(c) 


<0 




c)r"(c)=0,r"'(c) = 





c)r"(c)=0,r"'(c) = 







&rW(c) >0 




&r(4)(c) <0 




d)L<0 


d) r"(c)=0,r"'(c) = 


0, 


d)r"(c)=0,r"'(c) = 


0, 




r(4)(c)=0&L<0 




rW(c) =0&L<0 






2) c-3ot = 0: 










b) /'(c)=0&r"'(c) 


>0 







4 Further exact solutions 

The purpose of this section is not to discuss in detail the physical properties of the derived solutions; rather it is 
to show how to construct predefined solutions using the first five corollaries. The existence of a critical mass mc 
and the vaUdity the WEC at spatial infinity and that of the NEC will be the cases since for all solutions derived 
in this section and those discussed earlier, Eqs (|37] ) & (l46l) to ( |49l ). the function r(p) — p admits an expansion 
in powers of 1 /p at spatial infinity. We will use the inverse problem method by which we select a function 
r(p) and integrate Eqs (fTTT) & (fT4l) to obtain A(p) & 0(p), respectively, and derive V{p) from Eq ([T3l) . For the 
three solutions presented below, the expressions for A(p), V(p) & (^(p), if they are available, are given in the 
corresponding part of Appendix B. 

4.1 Black holes & naked singularities with V negative at spatial infinity 

Let r(p) be defined by 

r{p)=p--,c>0. (56) 

P 

Eqs ^T^, ([T3]l & (O lead to Eqs (IB. II ) for A, ^ &V. By corollaries 4 & 5, the solution is either a BH or a 
NS. Using Eq (l45l) with ai = —c^ & ^2 = we conclude that V is always negative at spatial infinity for positive 
masses {m> rric = 0) whatever the solution is. At the center c, A & V have opposite signs: A ~ (c — 3m) / (6p) 
&y ~-(c-3m)/(12cp2). 

Fig. |2] depicts a BH solution for c = 1 & m = 2/3. The horizon is at ph = 1.5806 (c < Ph < 3m by corollary 
2), the potential V crosses the axis at pi = 1.09995 and reaches its minimum value at p^in = 1.16804. The 
behavior of the solution is as stated in corollary 4. Fig.[3]depicts a NS solution for c = 1 & m = 1/3 the potential 
V of which is always negative. The behavior of the solution is as stated in corollary 5. 

4.2 Black holes & naked singularities with V having both signs at spatial infinity 

In the previous example (l56l ) it was not possible to derive a BH solution with V positive at spatial infinity or a 
NS with V positive near the center. We realize that by making another choice for r(p) 

r{p)=p-—^—,b>0,p>c = ^^^b>0, (57) 
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Figure 2: A BH solution with center c = 1 & mass m = 2/3 derived clioosing r(p) of tlie form l |56t . The horizon is at pj, = 1.5806 & 
the minimum of V is at = 1.16804. (^(c) = K. 




Figure 3: A NS solution with center c = 1 & mass m = 1/3 derived choosing r(p) of the form ( I56t . (^(c) = ;r & limp^c+ A(p) = 1/4. 

where b is any positive number and the center c = + 1)^/2. At spatial infinity we have r(p) = p — 
{b^ /A)/p — {b^ /A)/p^ + 0{l/p^). By corollary 6, the potential is positive as p ^ oo if m < = 3/4, re- 
gardless of the nature of the solution. A BH with c/3 < m < wic or aNS with m < min(c/3,mc) has its potential 
positive as p ^ oo. At the center c, on can show using Eqs (l40l ) & (|4T]) . with r(p) given by (|57] ). that A & V 
have always opposite signs for all m > & c > 0. 

Eqs (dB, dlS) & dJl lead to Eqs (IB.2I) to (IB. 41 ) for A, V & 0. By corollaries 4 & 5, the solution is either a 
BH or a NS. 

Taking m = 1/2 & c = (\/2 + 1)/2 (Z^ = 1), one obtains the BH solution depicted in Fig. |4] where the horizon 
is at ph = 1.31091 (c < ph < 3m by corollary 2). The potential crosses the axis at pi = 1.22754 & p2 = 2.46669 
and reaches its minimum and maximum values at Pmin = 1.2464 & p^ax = 2.81174, respectively. The solution 
behaves as stated in corollary 4. Taking m = l/3&c = l, one obtains the extreme NS depicted in Fig. [5] The 
function A is first concave up then concave down, which is obvious from the two plots of A in Fig. [5] The 
solution behaves as stated in corollary 5. 

4.3 A one-parameter particle-like solution 

In the two previous example (l56l) & (|57] ) it was not possible to derive a P-L solution. We realize that by making 
a third choice for r(p). We choose r(p) of the form (|37] ) with n = 4 (the simplest case). In order to fulfill 
the requirements of the theorem we take c = 3m. Since the Einstein-scalar theory is scale invariant, we restrict 
ourselves to the case c = 3m = 1. Once the expression of A(p) is obtained for c = 3m = 1, one generates the 
general expression of A for c = 3m / 1 by replacing p by p/c and m by m/c in the r.h.s of A(p). Hence, we 
start with 

4 1 fp-l)"^ 

r{p) = -{p-l)--^-^^,p>\. (58) 

Eq (ITTI ) leads to Eq (IB. 51 ) for A. The solution behaves as stated in the theorem with limp^f+ V{p) = — 31imp^c+ 
A"(p)/2 = —1.81735 (Eq (|35])). By corollary 6, positive masses (m > = —1/2) have negative potential at 
spatial infinity. 
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Figure 4: A BH solution with center c = (\/2+ l)/2 & mass m =1/2 derived choosing r(p) of the form {57}. The horizon is at 
Ph = 1.31091 and the minimum and maximum values of V are at Pn,in = 1.2464 & Pmax ~ 2.81174, respectively. (j){c) = 3;r/2. Two 
plots of V are shown for c < p < 1 .9 & 1 .9 < p < 5. 

The P-L solutions derived from (l37l) for « > 4 have similar behaviors with V negative everywhere. Plots of 
A & V for n = 5,6&1 and the same mass and center are shown in Fig. |7] As n oo^ A & V approach their 
Minkowskian limits. Such a behavior is desirable in theories where the effects of the scalar field are confined. 

5 Conclusion 

The selection criteria we have established in the first five corollaries and theorem have been used to help dis- 
cussing and constructing explicit solutions of particular interest. In section |2] we provided a unified formula for 
r(p ) (Eq (|37] )) which generates different P-L, NS & BH solutions (behavior of V and horizon location (for BHs) 
depend on n). In section [3] we have shown the existence of new kind of BHs (Eq ( |46l )) where the metric A is 
singular at the center only and the matter fields V & ^ aie regular everywhere. A parallel conclusion applies to 
the existence of NSs (Eq (|47])). 

We have shown that the NEC (any V) and the WEC (only V > 0) are always satisfied. Particularly, if r — p 
admits an expansion in terms of 1/p, then 1) there exists a sort of critical mass (if cJi 7^ 0) such that V at spatial 
infinity & {nic — m) have the same sign, regardless the nature of the solution. 2) The solutions to the field 
equations obey the WEC at spatial infinity (any V). If the derivatives of r at the center and L have specific signs, 
the WEC is satisfied there. 

We have constructed BHs, Fig.[2](resp. Fig.©, the potential function of which has one local minimum value 
(resp. one local minimum and one local maximum values). The two BHs behave differently at spatial infinity. 
We have also discussed NSs, Fig.[3](resp. Fig.|5]l, the potential function of which is negative everywhere (resp. 
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has both signs). We have constructed famihes (a family for each value of n > 4 in ( [37] )) of one-parameter P-L 
solutions, Figures [6] & |71 which have negative potential everywhere. Each family of P-L solution is different 
from that constructed in Ref [4]. 



Appendix A 

Appendix A is divided into two parts and its purpose is three-fold. In Part 1, where c = 3m and r(p) satisfies the 
conditions (19^ & we prove that a) hmp^f+ A'(p) is finite and b) if 15/'^ (c) -4/ (c)r'"(c) = (particularly 
if r"{c) = & r"'{c) = 0) then limp^c+ A"(p) is finite and provide a value for the latter. In Part 2, where c / 3m 
and r(p) satisfies the conditions (EUl & ^9^, we derive formulas (|40ll, (HB & (021) for A, V &(j). 



1) Part 1: c = 3m 

We introduce the variable x = p —c. While r(x+c) is supposed to be four times differentiable by hypotheses ( [T9l ) 
& (Unil, A'^\x + c) need be defined on a small half-open domain including the center ;c = 0. By Taylor's Theorem 
there exists a number < ^i{x) < x such that 

rix + c) = r'(c)x+ '^x^ + . (A.l) 
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V (n=5) V (n=6) V (n=7) 




Figure 7: P-L solutions with center c = 1 & mass ra = 1/3 derived choosing r(p) of the form l l37t with n = 5, 6&7 (left to right). 



Now, we introduce the continuous function 



r{x + c)/[/{c)x], ifx>0 
1, ifx = 



Since r{x + c) has only one zero: r(c) = 0, the function r{x) has no zero in the half-open domain x G [0,oo). 
The derivatives of t up to order 3 exist and are finite for x e [0,oo). Particulaiiy, the first derivative of T depends 
on r" and its second derivative depends on t'(0) = r"(c)/2r'(c) & t"(0) = r'"(c)/3r'(c). For allx > we 
have 

r(x + c) = r'(c)xT(x) . (A.2) 



Next, consider the continuous function a(x) defined by 
a(x) = 



xyr^{x + c) = \/[/\c)T\x)], ifx>0 
l/r'4(c), ifx = ■ 

The function a and its first three derivatives are finite and continuous on the domain x G [0,oo) and have no 
poles in it. By Taylor's Theorem there exists a number < ^2{x) <x such that 

a(x)- ' 2r-(c) 15r-^(c) - 4/(c)/-(c) <y"'{^2ix)) , .... 
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Notice how the derivative of order n — 1 of a depends on the derivative of order n of r. 

We introduce an arbitrary constant xq. <x <X() <oa and let Ci, C3, C2, 5 & A be the constants 

2 r° ^ d 2 n ^^^^ d r° ^ a r° d 

1 2/^(c) 15/^^(c)-4/(c)/"(c) S{x^) 

^ 2r'4(c)x2 r'5(c)xo 6/6(0) 6 ' 

15/^^ (c)- 4/ (c)/"(c) ^ _ 10/(c)/^(c)r(3)(c) - 15/^^(c) - /^(c)rW(c) 
~ 6/6(0) ' ~ 67^ ■ 

With c = 3m we rewrite ([TT]) using the new variable x and the constant Ci 

A=Cir2-2( T-j-^^ — ^d3')r^ {0<x<xo<^). (A.4) 



r4(3; + c) 

For X > (y > 0), 3;//(j + c) = a (3;) /y\ Thus using (IA.2I ) & (IA.3I ) in (IA.4I) leads to 

C, - 2C2 - ^^'l ^Hc).-rA') + ''''''' 



^3(0) ^ ^ 6r'6(c) 
with 5(x) = S{xo) + r a'"(^2 (3') ) dy . 

J.Vn 



Since a'" has no pole in the domain x G [0,°°), ^(x) is finite and continuous on it. 

Taking the derivative of (IA.5I ) with respect to x then the limit as x ^ 0+ we obtain Eq (l27l ). The non- 
vanishing contributions come from the second term in the r.h.s of (IA.5I ). which has a derivative 2tt'//^(c) ^ 
r"(c)//3(c) as X — > 0+, and the third one the derivative of which approaches —4r"{c)/r'^{c) as x — > 0+. We 
have thus proven that limp^c+ ^'{p) = ^'(c) = — 3r"(c) //^ (c) is finite whenever c = 3m. 

It is obvious from (IA.5I ) that a necessary condition for limp^f.+ A"(p) to be finite is 

I5r"\c)-4r'{c)r"'{c)=0. (A.6) 
We assume that (IA.6I ) holds. Now, A"(x + c) is provided by ^ 

„ (r2)"A-2 

A" = ^—L- , (A.7) 

with ruse (ICT): (r^)" = 2r'2(c)[T2 + 4xt't + x2t'2 +;(;2^^//] ^ ^^ ^^ 

All functions in the r.h.s of (IA.8I ) are continuous and finite for x e [0,oo) (t" depends on r'" which exists and 
is finite by hypothesis). Using (IA.2I ). (IA.5I ) & (IA.8I ) in (IA.7I ). one sees that the only non-finite contributions to 
limp^c+^"(p) may come from the second and third terms in (IA.5I ). These two contributions to A" are 

/ 3 8r"(c) _4^^ ^__2 4r"(c) , 4/'(c) 8r"(c) 



r(c) / r^[c)x r^[c)x r^[c)x 

Thus, we have proven that limp^c+ A"(p) is finite whenever condition (IA.6I ) is satisfied. 

As we have seen earlier, the potential V is finite at the center only if r"{c) = & r"'{c) = (Eqs (l30l) & 
(I35])). In this case Eq (IA.6I ) is satisfied. We postpone the derivation of a formula for A"(c) to Part 2. We will see 
that for c = 3m, r"{c) = & r"'{c) = , A"{c) reduces to L 

L = 2 hm f - 2/2 (c) r - ^ttVtI ' ^hi^h is just Eq m ■ (A.9) 
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2) Part 2: c^3m 

r{x + c) is supposed to be five times differentiable by hypotheses (1381 ) & ( [39l ). As before r^^^ (x + c) need be 
defined on a small half-open domain including the center x = 0. Replace dA.ll ) & ( IA.3I ). respectively, by 



r(x + c) =r'(c)jc + 
a(;c) 



2 

2/'(c) 



x^ + - 



(c) 3 , rW(c) 4 , 



""+^"+ 120 
aW(^2(x)) 



(a'^'^) (x = 0) depends on (c)). A (Eq (fTTI )) takes the form 

ydy 



A = -2r^(x + c) 



r^{y + c) 



'2{c — 3m)r (x + c) 
0(3;) d); 
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dy 



X , 



(A. 10) 
(A. 11) 



= [Ci-2(c-3m)C3]r2(x + c)-2r2(x + c) T-^^^ -2(c-3m)r2(x + c) /" 



o{y)dy 



Define the new continuous functions I\ (x) & I2 (x) by 



1 
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I,{x)=h{xo) + - ya^'\^2{y))dy, h{x) = h{xo) + - a^'\Uy))^y- 
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(A.12) 
(A. 13) 

(A. 14) 



Replacing o{y) in (IA.13I ) by the r.h.s of (lA.llI ) and using (IA.14I) . we arrive at 

1 lr"{c) 



A =[Ci-2{c-3m)C^y{x + c)-2r\x + c) 



h X 



+ ■ 



2r'4(c)x2 r'5(c)x 



+ 5 Inx + Ax — (xo x) 



2(c — 3m)r^{x + c) 



h{x) 



1 



+ 



/'(c) 



3r'4(c)x3 r'5(c)x2 x 



+ A Inx — (xo x) 



(A. 15) 



where (xq x) substitutes xq for x in the terms within the square parentheses. 

In order to obtain the behavior of A & V near the center (for small x), we replace r(x + c), /i(x) & h{x) 
in (IA.15I ) by the r.h.s of (lA.lOl ). /i(0) & /2(0), respectively, then extract the coefficients of 1/x, x, x^lnx, x^ 
and the independent term, we arrive first at ( |40l) . then using Q, we obtain (|4TI ). The derivation of (l42l ) is 
straightforward. The coefficient of x^ in (IA.15I ) reads 



A"{c) 



y-' r'2(c)xQ r'-'(c)xo 



1 . + i;:!(£L + 2/2(c)5lnxo 



+ 2r'"(c)Axo + 2/2 (c) [/i (xo) - h (0)] + 



r'"(c) 



15/'2(c) 

3/3 (c) ~ 4/4 (c) 



(A. 16) 



+ (c — 3m) 



2/2(c) 



-To 



g(y)dy 2 

3/2(c)xq /3(c)xj 



+ 



2/'(c) 2/2 (c) 5 



Xo 



2/2(c)Alnxo + 2r'2(c)[/2(xo) -/2(0)] +25r'(c)r"(c) + 



rW(c) /'3(c) 5/'(c)/"(c) 



18/3(c) 2/5(c) 9/4(c) 



A"(c) depends apparently on xq, however, using Eqs (lA.llI) & (IA.14I ) it is straightforward to obtain d[A"(c)]/dxo 
= 0. Hence, A"(c) is independent of xq and is finite since every function or constant in the r.h.s of (IA.16I ) is 
finite. In the limit xq ^ 0+, many terms drop from the r.h.s of (IA.16I) and the remaining terms provide the 
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desired expression for A"(c) 
A"{c) 



- lim 

.vo^O+ 



aiy)dy 



1 



+ (c — 3m) lim 



r 



■ + ■ 



4r"( 



-2r'^(c 



'■'^oa(y)d3; 2 



+ 2r'^(c)5ln;co 



(A. 17) 



+ 



^ ^ ^ ^ +2r'2(c)Alnxo 



Xo 



/"(c) 15r"2(c 



3r'3(c) 4r'4(c) 



+ (c — 3m) 



25r'(c)r"(c) + 



r(4)(c) r"3(c) 5r"(c)r'"(c) 



18r'3(c) 2r'5(c) 9r'4(c) 



where the limit in the r.h.s depends on the particular form of r(x + c) and cannot be obtained in closed form. If 
c - 3m = 0, r"(c) = & /"(c) = 0, Eq (IAT/I) reduces to (lAiOl) . 



Appendix B 

This section is divided into three parts. In Part 1 to Part 3, we show the expressions of A(p), V(p) & 0(p), if 
they are available, corresponding to subsections 14. 1 l to [4.3[ respectively. 

1) Part 1: 

The functions A(p), V(p) & 0(p) of subsections 14. 1 1 take the forms 

^ 3mp 11m c2 /c2-3mp\ ^ni{c^ - P^f\A^\ 
^ = l-^-^-^(,^2r^J+ ' = ^-2arccos(c/p), 



_ 15m ^ 5c^ - 6cV2 + 3mc2p3 3m(5c4 + Ac^p^ + 3p^) In 

^ ~ 8p3+8c2p^ 3p4(c2-p2)2 lec^p^ ' 

2) Part 2: 

The functions A(p), y(p) & 0(p) of subsections I4.2l take the forms 

^ _ (I8m-IlZ7)p 3Z^2_^42^^_2(Z,+ 18m)p [Wb - \%m)p^l{p) 
32^2 ^ 6(4p2-4Z^p-Z;2) + 32^2^3 

106Z^ + 36m + ^/2{\%m - 1 lb)^{p) 43b - 210m - 3^/2(11^? + 18m)£(p) 
^ 128^ ^ 384(p-Z^) 

264Z7m - 76^72 + 3^2^(1 IZ? - 18m)^(p) 
+ 3072(p-Z.)2 ' ^^-^^ 

25(llZ7-18m) 2[Z^2_^5^^_^2(ft-6m)p] 30ftm - 13^72 - 2(11^7 - 18m)p 

96b^{p-b) 3(4p2-4Z7p-Z72)2 + 4p2-4Z7p-Z72 

3(11^7- 18m)£(p) 10% -318m- 3^2(1 lZ7+18m)^(p) 

32V^^ ^ 192Z.(p-Z.)2 ^ ■ 

8% - 6m + 3V2{nb - lSm)i{p) 5[{264bm - 76b^ + 3^^^/2(11^^ - 18m)£(p)] 
384(p-Z^)3 3072(p-Z?)4 ' 

(B.4) 



+ 



^-4arctan((V2+l)>-^^+^)M 
2 r 2p + {V2-l)bJ 
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where Hp) = In f (1±^/^] = m I " 



{2-V2)p-b \ (V2+l)(p-c) 



3) Part 3: 



The function A(p) of subsections 14.31 is derived here. With our choice ( [58l) . the integrand in ([TTI ) is a rational 
function of the form -81piV[(p - 1)^(1 + ap)4(l + yp +j32p2)4]^ ^ ^ 2^/3 - 1, j8 = 2^3 + 1, y=_(22/3+2) 
& 7^ -4j82 < 0. With 77(p) = 1 + ap, a)(p) = 1 + yp + fi^p^ & = 4j32 - we obtain 



37rfcj32y 

107rj8^(2jj82-/7) n{2qa + hy) 2nvp^ ^nnj 6nlp^ a* (2 + 7P) ^ lOjjS^w' ^ /w' 



A? 2aA A3 A3 A^ A^co A^w 2A2w2 

A:(2 + 7p) 5jp^(o' Siyco' Sip^yco' 3l[5^(o' dkyo)' vco' jco' iyco' 
2A2w2 ^ 3A4w2 ~ 6A4w2 A^w ^ A^w ~ 2A4w ^ A2^ ^ 3A2w3 ~ 6j32A2w3 
, f 20P'{2jP^-ir) , 6PH2ll5^-kr) , 2(2v^2_ ) 2g« + /^7 ^ ^ ^«'Mo2,,, 
+ i A^ + + + ^^J^'^^^"UJ/^'^^' ^ ^ ^ 

where co' is the derivative of co with respect to p and a, e, f, h, i, j, k, I, p, jj., v & q are constants: 

a = - — , e = / = _?1±^, h = -— (75 - 58j8 + lljS^), / = — (-33 + 162^8 - 65p^), 

256 768 768 768 ^ ^ ^' 64 ^ h h j, 

j = ^ (^"^ - 82i8 +25/3^), ^ = ^ (12 -21/3 +7/3^), / = ^ (-62 + 71^8 - 19^8^), = A (3 -/3), 

M = (-123 + 128/3 - 54/3^), v = -i^ (167 - 120/3 +26^8^), q = ^ (-33 - 50^8 + 11^8^) . 

V is derived from ([T3] ) and ^ from (fT4l) . A computer-aided algebra leads to a sizable analytic expression for (p 
depending on the incomplete elliptic integral of the first kind and the complete elliptic integral of the third kind. 
Though the plot of (p shown in Fig.|6]has been drawn using the analytic expression of 0, a numerical evaluation 
of the integral in (fT4l) is as fast as the numerical evaluations of the two non-elementary elliptic functions of 
which depends. 
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